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ABSTRACT 

GRB 990123 was the first burst from which simultaneous optical, X-ray and gamma- 
ray emission was detected; its afterglow has been followed by an extensive set of radio, 
optical and X-ray observations. We have studied the gamma-ray burst itself as observed 
by the CGRO detectors. We find that gamma-ray fluxes are not correlated with the 
simultaneous optical observations, and the gamma-ray spectra cannot be extrapolated 
simply to the optical fluxes. The burst is well fit by the standard four-parameter GRB 
function, with the exception that excess emission compared to this function is observed 
below rsj 15 keV during some time intervals. The burst is characterized by the typical 
hard-to-soft and hardness-intensity correlation spectral evolution patterns. The energy 
of the peak of the vj u spectrum, E p , reaches an unusually high value during the first 
intensity spike, 1470 ±110 keV, and then falls to ~300 keV during the tail of the burst. 
The high-energy spectrum above ~ 1 MeV is consistent with a power law with a photon 
index of about -3. By fluence, GRB 990123 is brighter than all but 0.4% of the GRBs 
observed with BATSE, clearly placing it on the —3/2 power- law portion of the intensity 
distribution. However, the redshift measured for the afterglow is inconsistent with the 
Euclidean interpretation of the —3/2 power-law. Using the redshift value of > 1.61 and 
assuming isotropic emission, the gamma-ray fluence exceeds 10 54 ergs. 

Subject headings: gamma-rays: bursts 



1. INTRODUCTION the burst and also immediately after gamma-ray 

GRB 990123 was the first gamma-ray burst to emission was no longer detectable (Akerlof et al. 

be simultaneously detected in the optical band. 1999a,b). Observations of prompt X- and gamma- 

The Robotic Optical Transient Search Experi- ra y s were made h ? several spacecraft, including 
ment (ROTSE) detected optical emission during 
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all four instruments on the Compton Gamma- 
Ray Observatory, and the Wide Field Camera on 
Beppo-SAX (Feroci et al. 1999). The burst's 
afterglow was detected and monitored at X-ray 
(Heiseetal. 1999), optical (Odewahn et al. 1999) 
and radio (Frail et al. 1999) energies, resulting 
in the rapid localization of the burst (Piro et al. 
1999; Heise et al. 1999), the determination of the 
redshift (z > 1.61— Kelson et al. 1999; Hjorth 
et al. 1999), and the characterization of the af- 
terglow's spectrum and evolution (Galama et al. 
1999, Fruchter et al. 1999b, Kulkarni et al. 1999, 
Sari & Piran 1999b). This extensive set of obser- 
vations permits an analysis of the burst's relation 
to its afterglow. Here we present an analysis of 
the observations obtained with the instruments on 
the Compton Gamma-Ray Observatory (CGRO). 
Our goals are to: a) relate the X-ray and gamma- 
ray spectra to the optical flux; b) to characterize 
the burst, comparing it to other bursts, and c) to 
discuss the implications of this event. 

According to the currently favored class of the- 
oretical models, the observed burst and subse- 
quent afterglow radiation are synchrotron emis- 
sion by non-thermal electrons accelerated by 
strong shocks in relativistic outflows (Rees & 
Meszaros 1992; Galama et al. 1998; see Piran 
1999 for a review). In this model, the high-energy 
emission of the burst itself is radiated by "inter- 
nal shocks" which result from the collisions be- 
tween regions with different velocities (Lorentz 
factors) within the relativistic outflow, while the 
lower energy (radio, optical and X-ray) afterglow 
is thought to be emitted by an "external shock" 
where the outflow plows into the surrounding 
medium. Meszaros & Rees (1997) and Sari & 
Piran (1999a) pointed out that when the exter- 
nal shock first forms, a reverse shock propagates 
back into the relativistic outflow; the shocked re- 
gion behind the reverse shock is denser and some- 
what cooler than the region behind the external 
shock, and radiates in the optical band. The opti- 
cal emission from the reverse shock will be visible 
during or soon after the high-energy burst pro- 
duced by the internal shocks. 

Multi- wavelength observations during and im- 
mediately after the burst can test the current the- 
oretical model's predictions by detecting the tran- 
sitions between the emissions from the different 
regions. Soft X-ray tails, possibly the beginning of 
the X-ray afterglow, are sometimes observed after 
higher energy emission has ended (e.g., Yoshida 



et al. 1989; Sazonov et al. 1998). The X-ray 
emission at the end of GRB 970228 (Costa et al. 
1997) and GRB 980329 (in't Zand et al. 1998) 
is consistent with an extrapolation of the X-ray 
afterglow decay observed many hours after the 
burst, suggesting that the burst merges smoothly 
into the afterglow. GRAN AT/SIGMA observed 
a t-°- 7±0 03 power law decay of the 35-300 keV 
light curve for 1000 s following GRB 920723 (af- 
ter which the extrapolated flux would have been 
undetectable — Burenin et al. 1999); the power 
law begins immediately after the burst appears 
to end, when there is an abrupt change in the 
gamma-ray spectral index from to —1 (i.e., 
to Ne oc E~ x ). In prompt pointed observa- 
tions OSSE detected significant persistent emis- 
sion above 50 keV (> 300 keV in some cases) 
following the main burst in four of five events in 
a fluence selected sample (Matz et al. 1999). The 
emission was detectable for 10 2 -10 3 s with decays 
consistent with the power-law declines seen in the 
afterglows at lower energies. There was no evi- 
dence for a significant gap between the end of the 
burst and the beginning of the persistent emis- 
sion. BATSE observations of GRB 980923 show 
a power-law decay in time of the emission above 
25 keV between about 40 and 400 s after the be- 
ginning of the burst — this decay appears to be a 
higher-energy analog of the x-ray afterglows seen 
in other events (Giblin et al. 1999). On the other 
hand, the X-ray light curve of GRB 780506 ob- 
served by HEAO-1 shows a period of a few min- 
utes without detected flux after the burst, fol- 
lowed by renewed emission, suggesting a gap be- 
tween X-ray emission from the burst and from the 
afterglow (Connors & Hueter 1998). 

The ROTSE detection of optical flux simulta- 
neous with the gamma-ray emission and during 
the ~ 10 minutes following the burst (Akerlof et 
al. 1999a,b) can probe the issue of which physi- 
cal regions radiate when. These optical observa- 
tions consisted of three 5 s exposures beginning 
22.2, 47.4 and 72.7 s after the trigger of CGRO's 
Burst and Transient Source Experiment (BATSE) 
at 35216.121 s UT on 1999 January 23 (see Fig. 1), 
and three 75 s exposures beginning 281.4, 446.7 
and 612.0 s after the trigger. Gamma-ray emis- 
sion was detected for about 100 s, so the last 
three optical detections by ROTSE occurred after 
the apparent end of the burst phase of the GRB. 
ROTSE uses an unfiltered CCD, but an equiva- 
lent V band magnitude is reported; we will use 



M. S. Briggs et al. 



3 



these magnitudes without attempting any color 
corrections, etc. The optical flux increased from 
the first to the second exposure (m v = 11.82 to 
m v = 8.95), and then decayed from exposure to 
exposure (m v = 10.08, 13.22, 14.00 and 14.53). 

All four CGRO instruments detected this burst 
(however, the spark chamber of EGRET was not 
operating); in fact, ROTSE responded to a pre- 
liminary BATSE position. Here we describe the 
spectral evolution of the burst, and compare its 
spectrum with the optical flux measurements. In 
§2 we describe our analysis of the CGRO observa- 
tions of GRB 990123, while in §3 we discuss the 
implications of this analysis for burst models. 

2. DATA ANALYSIS 

2.1. Multi-instrument analysis 

The CGRO instruments observed GRB 990123 
at different energies and time scales. Here we 
present our procedures and some general results; 
more detailed discussion of the data of each in- 
strument appears in the following subsections. 

Figure 1 shows the light curve in different en- 
ergy bands accumulated by BATSE and COMP- 
TEL, the two instruments with the best combina- 
tion of temporal and spectral resolution for this 
purpose. The light curves show that the low- 
energy emission persists longer than high-energy 
emission — typical "hard-to-soft" spectral evolu- 
tion. Particularly striking is the virtual absence 
of high-energy emission 45-90 s after the burst 
trigger. In the highest energy band shown, 4 to 8 
MeV, emission is seen from the first spike but ap- 
pears to be absent or much reduced in the second, 
providing evidence that the first spike is harder 
than the second. The first spike also appears nar- 
rower in the 4 to 8 MeV band, even compared to 
the 2 to 4 MeV band. 

Fitting models to spectra provides the most ac- 
curate characterization of spectral evolution, al- 
though fits are possible only for spectra with a 
sufficiently large signal-to- noise ratio (S/N). For 
fitting the BATSE data we use the "GRB" func- 
tion (Band et al. 1993), which consists of a 
low-energy power law with an exponential cut- 
off, N E oc E a exp[— E{2 + a)/E p ], which merges 
smoothly with a high-energy power law, Ne oc 
E@. The "hardness" of this spectrum can be 
characterized by the energy E p of the peak of 
E 2 Ne oc vj v for the low-energy component (if 
P < -2, as is the case for GRB 990123, then E p is 
the peak energy for the entire spectrum). Most of 



the spectral curvature is found in the observations 
made with BATSE below about 1 MeV; the spec- 
tra above ~ 1 MeV from the other instruments 
are consistent with a simple power-law, which we 
identify as the high-energy power law E* of the 
GRB function. 

The model fitting is done using the standard 
forward-folding technique (see review by Briggs 
1996). Briefly, for each instrument we assume a 
photon model ('GRB' for BATSE, power law for 
the other instruments) and convolve that photon 
spectrum through a detector model to obtain a 
model count spectrum. The model count spec- 
trum is compared to the observed count spectrum 
with a goodness-of-fit statistic, and the photon 
model parameters are optimized so as to mini- 
mize the statistic. In most cases y 2 is use d as 
the goodness-of-fit statistic. This procedure is 
model dependent. Presenting a comparison of the 
data and model in count space is the best way to 
show the quality of the fit. However, it is diffi- 
cult to show results from all four instruments on 
one count-rate figure because of the widely dif- 
fering responses of the instruments. We therefore 
present photon spectra: the photon "datapoints" 
are calculated by scaling the observed count rate 
in a given channel by the ratio of the photon to 
count model rate for that channel; this ratio and 
therefore the photon datapoints are model depen- 
dent. 

An additional complication is that the differ- 
ing time boundaries of the spectral data from 
the CGRO instruments preclude forming a multi- 
instrument spectrum for exactly the same time 
interval. We have therefore selected time inter- 
vals which are as similar as possible and which 
include the two spikes of the burst (see Table 1). 
At high energies the flux of the spikes is much 
higher than the flux of non-spike intervals, con- 
sequently as long as each spectrum includes the 
spikes, differences in the time boundaries are rela- 
tively inconsequential, excepting that the spectra 
must be normalized to a common time interval 
encompassing the spikes. Because of the small 
differences expected from the differing time in- 
tervals and because of inter- detector calibration 
uncertainties, we have not attempted simultane- 
ously fitting the data from different instruments 
with the same photon model. 

The resulting photon spectra are shown in Fig- 
ure 2. The overlap between instruments is at the 
upper-end of the break region and in the high- 
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Fig. 1. — Top six graphs: Time history of GRB 990123 using data from CGRO detectors. The intervals over which 
the first three ROTSE observations were accumulated are indicated with shaded bands. Lower two graphs: The time 
evolution of E p , the energy of the peak of vf„, and of a, the low-energy photon index in the GRB function. These pa- 
rameters were obtained from GRB function fits to MER 16 channel spectra from BATSE LADs and 4. The horizontal 
bars on the points show the time over which the spectrum was fit, and the vertical bars show the la uncertainty in the 
values. 



M. S. Briggs et al. 5 



Table 1 

Multi-Instrument Spectral Analysis 







Inte 


gration Times 


Energy 


Range § 




Model Normalizationt 


instrument 


detector 


start^ 


endt 


length 


low 


high 




at 1 MeV 






(s) 


(s) 


(s) 


(MeV) 


(MeV) 




(7 cm" 2 s" 1 MeV" 1 ) 


BATSE 


LAD 


12.288 


45.056 


32.768 


0.033 


1.80 


-3.11 ±0.07 


1.93 ±0.02 


BATSE 


SD disc. 


12.288 


45.056 


32.768 


0.008 


0.013 






BATSE 


SD 1 disc. 


12.288 


45.056 


32.768 


0.015 


0.030 






BATSE 


SD 4 


12.736 


44.928 


32.192 


0.320 


25. 






OSSE 


Detector 2 


12.288 


45.056 


32.768 


1.0 


10. 


-2.82 ±0.16 


1.68 ±0.23 


COMPTEL 


BSA 


14.199 


46.398 


32.199 


1.0 


8.5 


-2.78 ±0.16 


1.49 ±0.10 


comptel 


telescope 


13.331 


46.099 


32.768 


0.75 


30. 


-3.33 ±0.26 


2.0 ±0.26 


egret 


TASC 


-0.057 


64.503 


64.560 


0.97 


250. 


-2.71 ±0.08 


1.09 ±0.07 


t Relative to BATSE trig 


ger time of 35216.121 s. 











s Energy range of the data used in the fit. 
"^Systematics dominate the statistical errors. 



energy power law regime. The overall agreement 
between the four instruments and the six detec- 
tor types is very good. Data from the BATSE 
detectors, LAD and SD 4, overlap from 0.32 to 
1.8 MeV and are in excellent agreement except for 
the highest energy LAD point at 1.1 to 1.8 MeV. 
The most discordant points are the highest en- 
ergy point from BATSE LAD and the two low- 
est energy points from the EGRET TASC. Be- 
cause of the very small statistical uncertainties 
on these three points, systematics in the calibra- 
tions probably dominate. The values of the high- 
energy spectral index (5 obtained by the instru- 
ments are largely consistent (Table 1). There are 
larger differences between the model normaliza- 
tions (Table 1). Differences in the normalizations 
at the level of ~ 10% are expected because of un- 
certainties in the effective areas. The assumption 
of a simple power-law for the OSSE, COMPTEL 
and EGRET fits, despite the presence of some 
curvature in the lowest energy data of these in- 
struments, also contributes to the discrepancies. 

2.2. BATSE 

The different BATSE datatypes provide vary- 
ing temporal and spectral resolution. Each of 
the eight BATSE modules contains a Large Area 
Detector (LAD), a relatively thin 2000 cm 2 Nal 
crystal, and a Spectroscopy Detector, a 12.7 cm 
diameter by 7.6 cm thick Nal crystal. The large 
circs, of the LADs permit high temporal resolution 
analyses, at limited spectral resolution, while the 
Spectroscopy Detector (SD) permit higher spec- 



tral resolution studies. The LADs are all operated 
with the same detector settings but the SDs are 
run at different gain settings and therefore cover a 
variety of energy ranges. Spectra are accumulated 
more frequently for the modules with higher count 
rates. Fortunately, the SD with the greatest count 
rate, SD 0, was in a high gain state, resulting in 
spectra from ~ 25 to 1750 keV ("SHERB" spec- 
tra), as well as a calibrated discriminator channel 
("DLSCSP1") between 8 and 13 keV. SD was at 
a burst angle of 35.2°. The second "rank" detec- 
tor, SD 4, was at a burst angle of 53.8°; with a 
low gain setting, SD 4 provided spectra between 
320 keV and 25 MeV. A discriminator channel be- 
tween 15 and 27 keV was obtained from the third 
rank detector SD 1 at an angle of 58.6°. 

The LADs provide data from about 33 to 
1800 keV in 16 channels. For Fig. 2, we use the 
CONT data from LAD 0, which provides 2.048 s 
resolution. The MER datatype from the LADs 
provide rates in the same 16 channels every 16 ms 
for the first 32.768 s of the burst, followed by 
64 ms resolution for an additional 131 s. For 
GRB 990123, the MER data is summed from 
LADs and 4, which had angles to the burst of 
27.5° and 46.0°, respectively. (For GRB 990123, 
high time resolution data from the energy channel 
from 230 to 320 keV is missing due to a telemetry 
gap; however, all channels are available at 2.048 s 
resolution via the CONT datatype). These MER 
rates show the burst's temporal morphology (see 
Figure 1), and are particularly useful for studying 
spectral evolution. 
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Fig. 2. — Deconvolved spectra from the CGRO detectors, shown both as photon flux Ne and in E 2 Ne = vf v units. 
The spectra have been rebinned into wider bins for clarity. Each spectra is calculated using the actual accumulation 
times (Table 1), except for the EGRET TASC spectrum, which uses a shorter time interval during which the emission 
was intense (see text) 
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Figure 1 (lower panels) shows the evolution of 
E p using fits to 16 channel MER spectra from 
LADs and 4 rebinned in time to provide S/N 
of at least 100. In order to improve the reliabil- 
ity of the fits, and because there is little evidence 
for temporal variations in (3, the GRB function 
was used with (3 fixed at —3.11. This value of 
(3 was obtained from the joint fit to the BATSE 
data shown in Fig. 2 and is consistent with the 
values obtained from the other instruments (see 
Table 1). As can be seen, E p increases by a large 
factor every time there is a spike in the light curve, 
as is typical of "hardness-intensity" spectral evo- 
lution. Additionally, the maximum E p is greater 
in the first spike than in the second, E p decreases 
more rapidly than the count rate and has an over- 
all decreasing trend, behaviors which are typical 
of "hard-to-soft" evolution (Ford et al. 1995). 
The small E p maximum for the second spike is 
consistent with the absence of that spike in the 4 
to 8 MeV lightcurve (Fig. 1). Two intervals dur- 
ing the first spike have E p values of 1470 ±110 
keV. Such values are exceptional — only 3 bursts 
of the 156 studied by Preece et al. (1999) have 
spectra with E p values above 1000 keV. 

To investigate the burst spectrum over the 
broadest energy range possible, we extend the 
LAD spectra by also fitting the SD data. The 
high-energy resolution SHERB data can be fit sat- 
isfactorily by the GRB function discussed above; 
the fits are consistent with the fits to other data 
types. SD 4 provides detections of burst flux 
to at least the 4.0 to 8.0 MeV band (Fig. 1). 
For the multi-instrument fit shown in Fig. 2, the 
BATSE data from LAD 0, SD 4 and SD discrimi- 
nators and 1 are jointly fit to a common photon 
model, allowing a small "float" in normalization 
between the detectors (for the discriminators, the 
relative normalizations are found using the high- 
resolution data from the same detectors). The 
relative normalization between LAD and SD 4 
is consistent with unity: 1.01 ± 0.01. The rel- 
ative normalization differences between the SDs 
are 10%, which is a typical value. The result of 
the fit to this 32 s interval is E p = 720 ± 10 keV, 
a = -0.60 ± 0.01 and (3 = -3.11 ± 0.07. At 
1 MeV, the fitted flux is 1.93 ±0.02 photons cm" 2 
s- 1 MeV" 1 . 

To test for the presence of an X-ray excess in 
GRB 990123 (as is seen in some other bursts — 
Preece et al. 1996), we fit the SHERB spec- 
trum and the discriminator channel accumulated 



by SD over various segments of the light curve, 
including each of the first three ROTSE observa- 
tions. We used the GRB function and fit all data 
points. For SD the discriminator channel covers 
8-13 keV, and the SHERB data 26-1760 keV. A 
significant soft excess above the fit of the GRB 
function is present in both of the major emission 
spikes (i.e., between 17.280 and 44.928 s after the 
BATSE trigger), including the first ROTSE ob- 
servation, but is not evident in tail of the burst 
(i.e., after 44.928 s after the trigger), including the 
second and third ROTSE observations. This soft 
excess is also not present during the weak emission 
before the first major emission spike (i.e., up to 
14.464 s after the trigger), and surprisingly, dur- 
ing a weak secondary peak on the declining edge 
of the first emission spike (29.696-34.304 s after 
the trigger). To verify the presence of the soft 
excess in the first ROTSE observation, we added 
the spectrum accumulated by SD 1 which was at 
a lower gain and therefore covered a higher en- 
ergy range. In particular, the SD 1 discriminator 
channel partially covers the gap between SD 0's 
discriminator channel and SHERB spectra. As 
shown by Figure 3, the joint fit is satisfactory, 
and the X-ray excess is weakly evident in the SD 1 
discriminator channel. 

Integrating spectral fits over the energy range 
20-2000 keV provides a measure of the intensity 
of the high-energy emission. During the three 
ROTSE observations the photon fluxes are 30.98, 
8.17 and 7.74 photons s" 1 cm -2 while the energy 
fluxes are 2.14 x 10" 5 , 1.68 x 10~ 6 and 1.82 x 10~ 6 
erg s -1 cm" 2 ; the ratios of these fluxes show that 
the average photon energy decreased from the first 
to the second and third observations. Integrat- 
ing the GRB function fit to the SHERB spec- 
trum accumulated over the entire burst gives a 
> 20 keV fluence of 3.0 x 10" 4 erg cm" 2 . The 
GRBM of Beppo-SAX reported a 40-700 keV flu- 
ence of 3.5(±0.4) x 10~ 4 erg cm" 2 (Feroci et al. 
1999). The > 20 keV fluence 5.09(±0.02) x 10" 4 
erg cm -2 reported by Kippen et al. (1999) was 
obtained using BATSE 4-channel discriminator 
data. This technique is less accurate, particu- 
larly for hard bursts like GRB 990123, because 
assumptions about the spectrum must be made 
in the deconvolution process. These diverse flu- 
ence values gives a sense of the consequences of 
differences in the inter-detector calibrations and 
analysis methods. 

The low-energy SD discriminator channel 



8 



CGRO Observations of GRB 990123 




Fig. 3. — Fit to SD count spectra for the time in- 
terval of the first ROTSE observation. The data points 
are from the SHERB data of SD and SD 4, and from 
the low-energy discriminators of SD and SD 1, as 
indicated. The model photon spectrum, the GRB func- 
tion, has been folded through the detector response, re- 
sulting in the model count spectra (histograms). Each 



detector, including each discriminator, has a separate 
count model and therefore a separate histogram. The 
histogram for each discriminator appears as a horizon- 

e tal line. Note that the lowest energy discriminator point 

I ; clearly lies above the GRB model (horizontal line) . The 

J — 1 — ' — spectral parameters of the model are E p = 1170 ± 30 

1000 10000 keV, a = -0.63 ±0.02 and /? = -4.1 ±0.4. 

Erergy (keV) 



(DISCSP datatype) extends our spectrum down 
to 8 keV for SD 0. The WFC on Beppo-SAX 
covers the energy range 1.5-26 keV. Feroci et al. 
(1999) report that the WFC flux peaked ~ 40 s 
after the high-energy peak, which would place 
the WFC peak between the second and third 
ROTSE observations. The WFC peak intensity of 
3.4 Crab corresponds to ~ 17 photons s _1 cm -2 
in their energy band. We fit the DISCSP and 
SHERB data for SD over ~ 9 s around the time 
of the WFC peak; only the low-energy part of 
the GRB function was necessary to fit the spec- 
trum. The DISCSP data point was ~40% higher 
(~ 1.5cr ) than the model fit, suggesting that there 
is at most a moderate soft excess. Integrating the 
model fit over the WFC energy range gives a flux 
of ~ 16 photons s _1 cm -2 , which is in acceptable 
agreement with the WFC flux, given the uncer- 
tain inter-detector calibration. 

2.3. OSSE 

The Oriented Scintillation Spectrometer Exper- 
iment (OSSE) consists of 4 collimated (3.8° x 
11.4°) Nal detectors covering the 50 keV to 
250 MeV energy range. The burst occurred ~ 20° 
and ~ 60° outside of the detector's narrow and 
broad aperture dimensions, respectively. How- 
ever, because the burst was intense and hard, 
radiation penetrated the shielding, permitting a 
> 1.0 MeV spectrum to be accumulated. In ad- 



dition, a high time resolution (16 ms) light curve 
was recorded above ~0.13 MeV by the Nal shields 
for the first 60 s of the burst. The time structure 
is complex with variations observed as short as 
~ 16 ms. Light curves can be formed in different 
energy bands using the photons which penetrated 
the shields into the detectors. These light curves 
show that most of the emission above a few hun- 
dred keV occurred during the two intensity spikes 
at the beginning of the burst; this is consistent 
light curves of BATSE and COMPTEL (Fig. 1). 
The spectrum appears to extend to the 3-10 MeV 
band. 

The spectrum of these photons can be analyzed, 
even though they did not originate in the de- 
tectors' field-of-view and were attenuated by the 
shields, because an appropriate response model 
was created and tested on solar flares at compara- 
ble off-axis angles. The spectrum from Detector 2 
suffered the least scattering and attenuation by 
material in the spacecraft. The background spec- 
trum uses data from 115 s before and after the 
burst. A power law model was fit to the data, re- 
sulting in a 1 MeV flux of 1.68±0.23 photons cm -2 
s- 1 MeV- 1 and a photon index (3 of -2.82 ±0.16; 
this fit was used to create the photon spectrum 
shown on Figure 2. Integrating this spectrum over 
the 32.768 s of most of the high-energy emission 
results in a fluence of (0.91 ±0.27) x 10 -4 erg cm~ 2 
from 1 to 10 MeV. 



M. S. Briggs et al. 



9 



2.4. COMPTEL 

COMPTEL detected with high significance 
the two main spikes of GRB 990123 in both 
its imaging telescope ("double scatter"; 0.75- 
30 MeV) and non-imaging burst-spectroscopy 
("single scatter"; 0.3-1.5 MeV and 0.6-10 MeV) 
modes. The burst-spectroscopy or BSA mode re- 
lies upon two Nal detectors in the bottom layer 
of the telescope (see Schonfelder et al. 1993 for 
instrument details). 

In imaging telescope mode, COMPTEL pro- 
vides detailed information on individual time- 
tagged photons with | ms time resolution. How- 
ever, for GRB 990123 these telescope data suf- 
fered from two limitations: 1) the telescope- 
mode effective area was at best a few cm 2 as 
GRB 990123 occurred nearly 60° from COMP- 
TEL's pointing direction; and 2) because of the 
event's high intensity, even at this low effective 
area a total of about 16% of this burst's live- 
time was lost due to telemetry limits (a maximum 
throughput of about 20 events per second). 

In contrast, the effective area in burst spec- 
troscopy mode was roughly two orders of mag- 
nitude greater than in telescope mode, and the 
deadtime was negligible. However, the spectral 
accumulation time was set at 4 seconds. There- 
fore the COMPTEL light-curves displayed in 
Fig. 1 are from the burst spectroscopy mode. The 
signal was negligible above 4 MeV. For both tele- 
scope and burst spectroscopy modes, the back- 
ground appeared stable so results depended little 
(a few %) on the choice of background integration 
time. 

Spectra from the two modes were handled dif- 
ferently. For the 0.75-30 MeV telescope data, one 
selects only events falling within a certain angle 
of the source position ("angular resolution mea- 
sure", or A.R.M.; Schonfelder et al. 1993), both 
reducing background and creating a nearly diag- 
onal response (e.g. Kippen et al. 1998). For 
the telescope spectrum displayed in Fig. 2, a 20° 
A.R.M. limit was used. The 32.768 s integration 
interval (see Table 1) was chosen both to cover all 
the significant gamma-ray emission and to allow 
the best live-time calculations. The background 
was taken from 131 seconds of data 15 orbits prior 
to the burst (see Kippen et al. 1998). These data 
were fit via the forward-folding technique assum- 
ing a simple power-law. Because of the few counts 
per bin and the background level, the goodness-of- 



fit statistic used was the background-marginalized 
Poisson statistic of Loredo (1992). The statistic 
indicates that this fit is good; however, the data 
are also consistent with a break at or below 1 
MeV (as indicated by the fit to the BATSE data) 
and also with a break at or above 6 MeV. The 
best-fit was 2.0 ± 0.26(£/lMeV)- 3 - 33±a26 pho- 
tons cm~ 2 s _1 MeV -1 , giving a total fluence of 
(1.0 ±0.3) x 10~ 4 ergs cm" 2 (0.75-30 MeV). Sys- 
tematic errors in the effective area and from dead- 
time are comparable in magnitude to the statisti- 
cal errors specified here. 

The single detector count spectra obtained in 
"spectroscopy mode" were processed as follows: 
the background was estimated from a spectrum 
of 140 seconds duration starting 202 seconds prior 
to the BATSE trigger. Eight high range (0.6-10.0 
MeV) detector spectra (4 s integration time each) 
covering a 32 s time interval (Table 1), were back- 
ground subtracted and summed. For these data, 
the forward-folding fitting was done using x 2 as 
the goodness-of-fit criterion. The location of the 
burst (zenith distance = 56.4°) results in an ef- 
fective detector area of 541 cm 2 , corresponding 
to 87% of the on-axis area. Assuming a single 
power law, the fit to the data from 1.0 to 8.5 
MeV (Fig. 2) gives best fit parameter values of 
normalization = (1.49 ± 0.10) photons cm~ 2 s _1 
MeV" 1 at 1 MeV and index (3 = -2.78 ± 0.16. 
The fluence (1.0-8.5 MeV, 32 s) is 7.93 x 10" 5 erg 
cm~ 2 We note a clear break of the data below 1 
MeV where the spectrum becomes flatter. Pre- 
liminary analysis of the low range (0.3-1.5 MeV) 
spectroscopy data (not shown in Fig. 2) covering 
the same 32 s time interval also indicates a spec- 
tral break at around 0.8 to 1 MeV. 

2.5. EGRET 

The Energetic Gamma-Ray Experiment Tele- 
scope (EGRET) has two methods of detecting 
gamma-ray bursts. The EGRET spark chambers 
are the primary detector for the telescope and are 
sensitive to gamma rays above 30 MeV. The spark 
chamber detector is turned on in response to an 
onboard trigger provided by BATSE when a burst 
is detected within 40 degrees of EGRET's princi- 
pal axis. However, GRB 990123 was too far off 
axis (56.4°), so no information is available from 
the spark chambers for this burst. 

EGRET also detects gamma-ray bursts with 
the Total Absorption Shower Counter (TASC), a 
77 x 77 x 20 cm Nal(Tl) detector. While acting 
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principally as the calorimeter for the spark cham- 
bers, it is also an independently triggered detec- 
tor in the energy range from 1 to 200 MeV. The 
TASC is unshielded from charged particles and 
has a high background rate, which is determined 
from the time intervals just preceding and follow- 
ing the burst. The response of the TASC varies 
strongly with incident angle, because of the in- 
tervening spacecraft material. An EGS-4 Monte 
Carlo code is used with the CGRO mass model 
to determine the response of the detector 
function of energy and incident photon direction. 
In normal mode a spectrum is accumulated every 
32.768 seconds. Shorter time intervals are also ac- 
cumulated following a BATSE burst trigger; how- 
ever, these short intervals were prior to the peak 
of this burst since BATSE triggered 25 seconds 
prior to the peak of the emission. 

The burst was detected in the EGRET TASC in 
two successive 32.768 second intervals. The first 
is from -0.057 to 32.711 s relative to the BATSE 
trigger, the second from 32.735 to 65.503 s. The 
break between the two records occurs roughly 
between the two major spikes. The differen- 
tial photon spectrum is well fit by a power law, 
F(E) = Zc(i?/lMeV) /3 , in each separate time in- 
terval as well as when the intervals are combined. 
The combined 65.5 second interval (Fig. 2) is fit 
by (3 = -2.71 ± 0.08 and k = 1.09 ± 0.07 photons 
cm" 2 s" 1 MeV" 1 . The normalization k and the 
points in Fig. 2 are calculated assuming that the 
excess flux above background was emitted over 
32.7 s, the time interval of the high energy emis- 
sion from the burst as seen in the light curves of 
BATSE and COMPTEL shown in Fig. 1. 

The separate spectral fits of the two major 
spikes show no evidence of spectral evolution. The 
differential photon spectral index, /3, for the first 
spectrum is —2.56 ± 0.11, and for the second is 
-2.78 ± 0.15. The normalization k is 0.97 ± 0.08 
for the first spectrum and 1.20 ± 0.13 for the sec- 
ond spectrum, where the time intervals for the 
emissions of the first and second spectra are as- 
sumed to be 20.4 s and 12.3 s, respectively. 

3. DISCUSSION 

3.1. Optical 

The detection of simultaneous optical emis- 
sion with ROTSE is the primary distinction of 
GRB 990123, and the obvious question is the 
relationship between the optical and gamma-ray 
fluxes. The optical and gamma-ray light curves 



are clearly very different: the optical flux rose by 
an order of magnitude from the first to second ob- 
servation, and then fell by a factor of 2.5 to the 
third one. On the other hand, the first optical ob- 
servation occurred during the first hard gamma- 
ray spike, while the second and third observations 
occurred during the soft tail after the two main 
gamma-ray spikes. Indeed, the 20-2000 keV en- 
ergy flux, calculated by integrating over fits to the 
spectrum, is an order of magnitude higher during 
the first optical observation than during the sub- 
sequent two, which had comparable fluxes. Thus 
the optical and high-energy fluxes are not directly 
proportional. 

But perhaps the optical and gamma-ray emis- 
sion are part of one component, and the rela- 
tive fluxes change as the spectrum evolves? Fig- 
ure 4 shows the observed optical and gamma-ray 
spectra for the three observations. The curves 
show the photon models obtained from fitting the 
BATSE data. Note the extent of the X-ray ex- 
cess on Figures 2, 3 and 4. As can be seen, the 
gamma-ray spectra do not extrapolate down to 
the optical band without an unseen break or bend. 
We have included the maximum 3a variations in 
the low-energy gamma-ray spectrum permitted 
by the data; these extrapolations underestimate 
the optical flux by an order of magnitude or more. 
Perhaps the optical emission and the X-ray excess 
observed during the first ROTSE observation are 
part of the same component? We extrapolated a 
power law with index a ox = —0.243 connecting 
the optical and 10 keV X-ray fluxes during this 
first observation. Figure 4 shows that a similar 
power law over predicts the X-ray flux by more 
than an order of magnitude during the second 
and third optical observations. We conclude that 
the optical emission is in excess of the low-energy 
extrapolation of the burst spectrum, which sug- 
gests that the optical and high-energy emission 
originate from different shocks in the relativistic 
outflow. 

3.2. X-Ray Emission 

Also interesting is the observation that the X- 
ray excess and the gamma-ray emission are ap- 
parently not related. That is, the excess occurs 
during the two main spikes in the gamma-ray time 
history, but not in the time interval afterwards, as 
well as being absent during a small flaring episode 
on the falling portion of the first main spike. The 
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X-ray time history in the 8-13 keV band shows 
interesting behavior during the tail end of the 
burst: rather than the flat shoulder present in the 
gamma-rays, there is a somewhat steady fall-off in 
flux from T+50 s to T+130 s. At the same time, 
the low-energy power-law index a makes a tran- 
sition from ~ —0.5 to ~ —1.0 (see Fig. 1). Thus, 
the portion of the burst associated with the X-ray 
excess is also associated with a hard low-energy 
power-law index; when the low-energy index soft- 
ens, the excess disappears. A possible reason is 
that the change in slope a causes the X-ray ex- 
tension of the gamma-ray spectrum to mask the 
low-energy excess. However, the flux predicted 
by the GRB function fit to the data of the second 
ROTSE interval is inadequate to mask the X-ray 
excess observed during the first ROTSE interval. 
This suggests that the X-ray component has dis- 
appeared by the second half of the burst — making 
GRB 990123 the first burst in which BATSE has 
reported evolution in this component. In that 
case, it is interesting that the X-ray excess is not 
correlated with the optical emission, which peaks 
during the second half of the burst. 

There are several burst models that support 
the existence of excess X-rays simultaneously with 
prompt burst emission. The most successful of 
these, the Compton attenuation model (CAM— 
Brainerd 1994), predates the discovery of the ex- 
cess. In this model, a beam of photons is scat- 
tered out of the line of sight by an optically-thick 
medium via Compton scattering in the Klein- 
Nishina limit, which is less efficient at higher pho- 
ton energies than Thompson scattering. It is this 



Fig. 4. — Optical flux and gamma-ray spectra for the 
time intervals of the 3 ROTSE observations. The curves 
for the first (third) time interval have been shifted up 
(down) by a factor of 1000. The data points "*" on the 
left are the ROTSE observations, while the solid curves 
on the right are the best fit models to the BATSE SD 
spectra. The dashed curves on the low-energy end of 
the gamma-ray spectra show the spectra for 3<r varia- 
tions in a, the GRB function's low-energy spectral in- 
dex. The dot-dashed curve is the extrapolation of the 
optical flux to 10 keV using a ox , the power law index 
connecting the observed optical and 10 keV X-ray fluxes 
for the first ROTSE time interval. 



change in efficiency that modifies an assumed in- 
put power-law spectrum into the observed shape, 
with a universal break energy in the rest frame 
determined by the physics. The low-energy up- 
turn is determined by the density of the inter- 
vening scattering material. In the simplest form 
of the CAM, there is no allowance for changes in 
the optical depth as a function of time, yet the 
change in amplitude of the X-ray upturn in the 
present burst seems to require this. When we fit 
the BATSE spectral data for the first ROTSE in- 
terval with the CAM, a x 2 value of 511 for 400 
degrees of freedom is obtained. Values of reduced 
X 2 above unity are typical of very bright bursts 
like GRB 990123 because systematic deviations 
dominate the statistical fluctuations. The best-fit 
value for the red-shift is 1.19 ±0.06; while x 2 in- 
creases by 27 to 538 when the red-shift is required 
to be 1.61. When the GRB model is fit to all of 
the BATSE data, including the SD discriminator 
points with their evidence for an X-ray excess, 
the same x 2 value of 538 is obtained. If the SD 
discriminator data are ignored, the GRB model 
has a substantially better x 2 value than the CAM 
model, indicating the power of these data for test- 
ing models with different X-ray predictions. 

Other burst models have various levels of com- 
patibility with an excess of X-rays, but few re- 
quire it. Usually, a separate, lower-energy spec- 
tral component would indicate that Compton up- 
scattering is producing the observed gamma-rays 
from a seed population of X-rays. One would 
expect the two components to be highly corre- 
lated in time, which is not observed here. The 
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input spectral shape is modified by the distri- 
bution of energetic particles that are doing the 
up-scattering, so one can infer a power-law index 
for the particles from the photon continuum spec- 
tral index. Finally, it is also not clear if there is 
any connection between the X-ray flux and the 
burst afterglow, presumably created in this case 
by a shock formed when relativistically-expanding 
material from the burst meets material external 
to the system, such as interstellar material, or 
the remnant of a wind from the burst's progen- 
itor (Sari & Piran 1999a). Although the compo- 
nent of the aftershock related to the optical flash 
(the reverse shock) is not expected to produce any 
detectable X-ray component, the forward shock 
is. The observation that the hard X-rays (top 
lightcurve of Fig. 1), have a maximum after the 
occurrence of the two main spikes in the gamma- 
rays, may indicate the peak of the forward shock 
in this band. 

3.3. Afterglow 

Park et al. (1997a) introduced three dimen- 
sionless statistics to compare optical detections or 
upper limits with the observed gamma-ray emis- 
sion. Behind these statistics are implicit assump- 
tions about the relationship between the emis- 
sion in these two bands. The first two statistics 
can be calculated only if the high-energy emis- 
sion is still detectable at the time of the opti- 
cal emission. Comparing the optical and gamma- 
ray flux densities f v (in units of erg s" 1 cm -2 
Hz" 1 ), is Ri = fu{vo P t)/fv{v-y), where v opt and 
v 1 are fiducial optical and gamma-ray frequencies 
(here the V band's 5.6 x 10 14 Hz and 100 KeV 
= 2.4 x 10 19 Hz). Comparing the optical flux to 
the gamma-ray flux density / 7 (erg s" 1 cm" 2 ) is 
Ri = fv(vapt)vopt/f"e Finally, R 3 = /^(^pt)/^ 
compares the optical flux density to the gamma- 
ray energy fluence 7 (erg cm" 2 ) from the be- 
ginning of the burst until the time of the op- 
tical observation. i? 3 is relevant and calculable 
whether or not the high-energy emission is de- 
tectable at the time of the optical observations; 
however the time between the high-energy emis- 
sion and the optical observation is physically in- 
teresting. While these statistics do depend on 
the fiducial optical and gamma-ray frequencies 
which are chosen because of detector capabilities, 
they do not depend on instrumental quantities 
such as the length of the optical exposure. Ta- 
ble 2 shows the values of these 3 statistics for 



the 6 ROTSE detections of GRB 990123. Also 
shown are typical upper limits obtained with the 
Gamma-Ray Optical Counterpart Search Exper- 
iment (GROCSE), ROTSE's less-sensitive pre- 
decessor, and the upper limits for GRB 970223 
(Park et al. 1997b) obtained with the Livermore 
Optical Transient Imaging System (LOTIS), an 
instrument very similar in design and operation 
to ROTSE (LOTIS did not attempt to observe 
GRB 990123 because of bad weather— H.-S. Park, 
personal communication 1999). As can be seen, 
the ROTSE detections are below the GROCSE 
upper limits, but are comparable to the upper 
limits for GRB 970223. Only with more detec- 
tions and upper limits from ROTSE, LOTIS and 
similar instruments will we determine whether 
GRB 990123's optical emission was typical. 

The i?3 statistic can also compare the optical 
afterglows to the gamma-ray emission; Table 3 
shows i?3 for R band afterglow emission after 
1 day for a number of recent bursts. These val- 
ues of i?3 are usually smaller than the detections 
and upper limits during and immediately after 
the burst because it is possible to observe the 
optical transient with large telescopes once the 
position is known accurately. We would like to 
use i?3 values calculated at the same time and 
optical frequency in the rest frame of the burst. 
Assuming that the afterglow flux has power law 
spectrum and temporal decay, f v oc u e t s , then 

^burst = ^observer t 1 + z Y~ e - For most af- 
terglows e ~ 5 ~ — 1 (e.g., Galama et al. 1997, 
1998, 1999; Bloom et al. 1998), and by afterglow 
theory 5 — e = (5 + l)/3 for adiabatic cooling; 
thus R 3 b^-gt has at most a weak dependence on 
redshift. As can be seen, there is a wide range 
of values, which confirms that afterglows are not 
all alike. Indeed, the i? 3 value and the isotropic 
gamma-ray energy are inversely correlated, which 
results from a broader energy distribution than 
the distribution of R band specific luminosities in 
the bursts' frame (i.e., the optical flux converted 
into a luminosity emitted by the afterglow). This 
suggests that the efficiency with which the to- 
tal energy is converted into emission varies with 
wavelength band. For example, in the standard 
model, the energy release by internal shocks de- 
pends on details of the relativistic flow, which can 
vary from burst to burst, while the kinetic energy 
of the flow is expended at the external shock; thus 
the distribution of afterglow energies may be nar- 
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Table 2 

Optical-Gamma-Ray Emission Comparison 



Obs. a 


t b 


M^opt) C 


/fiC-ElOOkeV 


) d 


/7° 




W 






R2 h 




Rs 1 




GRB 990123-1 


22.2 


6.79 x 10~ 25 


1.07 X 10" 


-8 


2.14 x 10 


-s 


7.78 x 10" 


-5 


1.53 x 10 1 


1.78 x 10 




8.73 x 10" 


■ 21 


GRB 990123-2 


47.4 


9.55 x 10" 24 


3.93 x 10" 


-9 


1.68 x 10 


-6 


3.05 x 10" 


-4 


5.87 x 10 2 


3.19 x 10 


-3 


3.13 x 10" 


■20 


GRB 990123-3 


72.7 


3.37 x 10- 24 


4.13 x 10" 


-9 


1.82 x 10 


-6 


3.39 x 10" 


-4 


1.97 x 10 2 


1.04 x 10 


-3 


9.95 x 10" 


-21 


GRB 990123-4 


281.4 


1.87 x 10~ 25 










2.65 x 10" 


-4 








7.06 x 10" 


-22 


GRB 990123-5 


446.7 


9.12 x 10" 26 










2.65 x 10" 


-4 








3.44 x 10" 


-22 


GRB 990123-6 


612.0 


5.60 x 10" 26 










2.65 x 10" 


-4 








2.11 x 10" 


-22 


GROCSEJ 


~ 25 
















< 1.5 x 10 4 


< 4 x 10 


-2 


< 3 x 10" 


-18 


GRB 970223 k 


11 


< 1.4 x 10" 24 


1.4 x 10" 


-9 


6.09 x 10 


-7 


4.8 x 10" 


-5 


< 2.4 x 10 2 


< 1.3 x 10 


-3 


< 2.9 x 10" 


■20 



a Optical observation. 



Time in seconds from BATSE trigger to beginning of optical observation. 
C V band optical flux density (erg s" 1 cm" 2 Hz -1 ). 
d 100 keV flux density (erg s" 1 cm" 2 keV" 1 ). 
e 20-2000 keV gamma-ray energy flux (erg s" 1 cm" 2 ). 

f 20-2000 keV energy fluence from the beginning of the burst to the middle of the optical observation (erg cm" 2 ). 
g Ratio of optical to gamma-ray flux densities, R\ = 2.42 x 10 17 /, y (i' pt)//£;(^iookev) (the constant is necessary for 
units conversion). 

h Ratio of optical flux density to gamma-ray energy flux, R 2 = fu{v -pt)v ov tl '/ 7 - 
'Ratio of optical flux density to gamma-ray energy fluence, R% — fv(yopt)/4>T 
j Typical GROCSE upper limits from Park et al. (1997a). 
k From Park et al. (1997b). 



Table 3 
Burst Energetics 



GRB 


z 


4><y (> 20 keV) a 


E b 


m R (1 day) c 


fA"R) d 


As 




Rcf. f 


970228 




2.0 x 10" 6 




21.1±0.15 


1.12xl0" 5 


5.6x10" 


23 


1,2 


970508 


0.835 


3.7xl0- 6 


5.3xl0 51 


20.85±0.05 


1.32xl0- 5 


3.8x10" 


23 


3,4,5 


971214 


3.418 


l.lxlO" 5 


2.5xl0 53 


23.0±0.22 


1.94X10" 6 


1.8x10" 


■24 


6,7 


980326 




2.2 x 10" 6 




23.39±0.12 


1.36X10" 6 


6.2x10" 


■24 


8 


980329 


~5 


5.0xl0" 5 


2.4 xlO 54 


23.9±0.2 


8.48X10" 7 


1.7x10" 


■25 


9,10 


980519 




1.1 x 10" 7 




21.28±0.13 


9.47X10" 6 


8.6x10" 


■22 


11,12 


980613 


1.096 


1.7xl0" 6 


4.2xl0 51 


23.2±0.5 


1.61X10" 6 


9.5x10" 


■24 


13,14 


980703 


0.966 


4.6xl0" 5 


8.9xl0 52 


21.18±0.10 


1.04X10" 5 


2.3x10" 


■24 


15,16 


990123 


1.61 


3.0xl0" 4 


1.6xl0 54 


20.4±0.3 


2.13X10" 5 


7.1x10" 


■25 


this paper 


a Encrgy fluence, in erg cm 2 . 

b Total energy in ergs if radiated isotropically, 


assuming H 


= 70 km s" 


1 Mpc" 1 


) ^0 


= 0.3, 



and A = 0. 

C R band magnitude after 1 day. 

d Equivalent R band flux density, in Jy. 

c Ratio of optical flux density to gamma-ray energy fluence, i? 3 = fv{ v R)/4>-y 
'References: 1. Galama et al. (1997); 2. Fruchter et al. (1999a); 3. Galama et al. (1998); 
4. Metzger et al. (1997); 5. Bloom et al. (1998); 6. Diercks et al. (1997); 7. Kulkarni et 
al. (1998); 8. Groot et al. (1998); 9. Palazzi et al. (1998); 10. Fruchter (1999); 11. Vrba et 
al. (1998); 12. Djorgovski et al. (1998b); 13. Djorgovski et al. (1998c); 14. Djorgovski et 
al. (1999); 15. Vreeswijk et al. (1999); and 16. Djorgovski et al. (1998a). 
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rower than that of the burst itself. 

3.4. Intensity 

GRB 990123's fluence 7 = 3.0 x 10~ 4 erg cm- 2 
places it firmly on the N(> 7 ) oc 0~ 3//2 por- 
tion of the fluence distribution — only 0.4% of the 
bursts observed by BATSE have had higher flu- 
ences. By photon peak flux, only 2.5% of BATSE 
bursts have been more intense. Whatever the in- 
tensity measure, peak photon flux P or fluence 
7 , the high end of BATSE's cumulative inten- 
sity distribution is a power law with an index of 
—3/2. In the simplest cosmological model stan- 
dard candle bursts occur at a constant rate per co- 
moving volume; bursts in nearby Euclidean space 
produce the 0~ 3 ^ 2 portion of the cumulative peak 
flux distribution while the bend in this distribu- 
tion results from bursts at sufficiently high red- 
shift (z > 0.3) that spacetime deviates from Eu- 
clidean. Thus we would expect GRB 990123 to 
occur at low redshift. However the absorption 
lines in the optical transient's spectrum show that 
z > 1.61 for GRB 990123 (Kelson et al. 1999, 
Hjorth et al. 1999). GRB 990123 demonstrates 
that the 0~ 3//2 portion of the peak flux distribu- 
tion does not originate only in nearby Euclidean 
space. Consequently, the burst rate per comoving 
volume must vary with redshift, with the burst 
rate and cosmological geometry "conspiring" to 
produce the Euclidean-like 0~ 3 ^ 2 dependence. In- 
deed it has been suggested that the burst rate 
is proportional to the star formation rate (e.g., 
Wijers et al. 1998; Totani 1997), which has de- 
creased precipitously since the epoch correspond- 
ing to z ~ 1.5. 

The redshift implies that GRB 9901 23 's energy 
release was at least 1.6 x 10 54 erg (H = 70 km s" 1 
Mpc -1 , Q = 0.3 and A = 0), if the burst radiates 
isotropically, comparable to the energy release in 
GRB 980329, if z ~ 5 for this burst (Fruchter 
1999). Energy releases are given as isotropic 
equivalents, even though there are indications of 
beaming (Fruchter et al. 1999b, Kulkarni et al. 
1999, but see Meszaros & Rees 1999) because the 
beaming angles are unknown. See Table 3 for a 
list of the energies of the bursts with redshifts. 
Both GRB 971214 (z = 3.418— Kulkarni et al. 
1998) and GRB 980703 (z = 0.966— Djorgovski 
et al. 1998a) had energies of ~ 10 53 erg while 
GRB 970508 (z = 0.835— Metzger et al. 1997; 
Bloom et al. 1998) and GRB 980613 (z = 1.096- 
Djorgovski et al. 1999) had energies of ~ 6 x 10 51 



erg. Thus we see that the standard candle as- 
sumption behind the simplest cosmological model 
is violated, and the intrinsic energy distribution 
is very broad. 

Thus two basic premises of the simplest cosmo- 
logical model are violated. The burst rate and 
the energy release distribution must now be de- 
termined empirically, with the peak flux distribu- 
tion providing a constraint on these distributions. 
Note that it may not be possible to separate the 
burst rate into independent functions of redshift 
and energy. 

4. SUMMARY 

The high-energy emission of GRB 990123 was 
fairly typical of gamma-ray bursts. The spectrum 
above 25 keV can be fit satisfactorily by the sim- 
ple four-parameter GRB function. The parame- 
ter E p , the energy at which vf v peaks, reached 
the unusually high value 1470 ±110 keV. The 
typical spectral evolution patterns of hardness- 
intensity correlation and hard-to-soft evolution 
are observed. An X-ray excess is present below 
~ 15 keV during the early part of the burst, 
particularly during the first ROTSE observation. 
The burst's peak flux and energy fluence place 
the burst in the top 2.5% and 0.4% of all bursts 
observed by BATSE; these values also are clearly 
in the portion of the intensity distribution which 
appears consistent with the —3/2 slope power law 
expected for sources uniformly distributed in Eu- 
clidean space. However, GRB 990123's redshift of 
z > 1.61 shows that the burst originated far be- 
yond nearby Euclidean space. This redshift and 
the fluence imply an energy release of > 1.6 x 10 54 
ergs (H = 70 km s" 1 Mpc" 1 , = 0.3 and 
A = 0), if radiated isotropically. GRB 980329, 
if at z ~ 5 (Fruchter 1999), has a similar en- 
ergy release; the determination of the redshift of 
GRB 990123 via spectroscopic lines in the after- 
glow leaves no doubt that isotropic-equivalent en- 
ergy releases above 10 54 ergs occur. 

The gamma-ray and ROTSE optical fluxes are 
not correlated, and the gamma-ray spectrum does 
not extrapolate down to the optical observations. 
Thus there is no indication that the emissions in 
these two bands are from the same component. 
Indeed, the standard burst theory attributes the 
high-energy emissions to "internal shocks" within 
a relativistic outflow, and the simultaneous opti- 
cal emission to a reverse shock that forms when 
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the relativistic outflow plows into the surrounding 
medium. 

Various measures that relate the optical and 
gamma-ray emissions show that the ROTSE de- 
tection of optical emission from GRB 990123 is 
consistent with the upper limits on simultaneous 
optical emission from previous bursts. Therefore, 
only future optical observations by ROTSE, LO- 
TIS and similar robotic systems will determine 
whether GRB 990123 was typical. 
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